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ABSTRACT 

SETTLEMENT PREDICTIONS BY STANDARDIZED METHODS 

by 

WILLIAM J. ANDERSON, JR. 



Submitted to the Department of Civil Engineering 

on August 13, 1973 

in partial fulfillment of the requirements 
for the degree of 

Master of Science in Civil Engineering 



The purpose of this thesis is to investigate the 
feasibility of using standardized methods of predicting total 
settlement and time rate of settlement and to evaluate the 
role of engineering judgement in these predictions. 

Several standardized methods have been presented 
and standardized procedures for selecting parameters have 
been recommended. Settlement was predicted for seven field 
cases using the various methods presented. The resulting 
predictions were then compared with actual measured perfor- 
mance. A valid comparison between observed and predicted 
long term settlement was found impossible unless the average 
degree of consolidation was known. A comparison of initial 
settlements was feasible. Results indicated that standardized 
methods are feasible, but that standardized procedures for 
selecting parameters are not. Engineering judgement is re- 
quired in the selection of the appropriate parameters to 
match the specified method. The results also showed the 
need for sufficient representative data and the need for a 
reasonable method of comparison between observed and pre- 
dicted performance. 
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CHAPTER I 
INTRODUCTION 



1.1 PREDICTIONS 

It is the job of the Civil Engineer to make 
predictions. A good prediction can not only save a client 
large sums of money, but may also help to prevent injury or 
loss of life due to failure. An engineer may be required' to 
predict the: 

1. Completion Date; 

2. Construction Cost; 

3. Factor of Safety; 

4. Environmental Consequences; 

5. Performance of the Structure. 

As part of the performance predictions, a geotechnical engineer 
may have to forecast the magnitude, direction and rate of de- 
formation of the structure on its foundation; the stability 
of the structure on its foundation; and the loads, both lateral 
and vertical, applied to structures resting on or within the 
soil (Lambe, 1973a). In order to make a prediction in a soil- 
rock system, the engineer uses the following steps: 

1. Determines the field situation and 
then simplifies; 
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2. Determines the mechanism (s) ; 

3. Selects a method and parameters; 

4. Manipulates the method and parameters to 
obtain a prediction; 

5. Portrays the prediction (Lambe, 1973a) 

Although all phases of the prediction process are important, 
probably the most crucial are steps 3 and 4. The relationship 
which exists between the method and parameters is vital in 
making a good prediction. 

The selection and manipulation of method and para- 
meters in a soil-rock system is complicated by many factors. 
Some of these are listed below: 

1. One cannot see underground; 

2. Soil-rock strata may vary widely 
both horizontally and vertically; 

3. Soil is not an isotropic and elastic 
medium; 

4. Soil properties are dependent on 
stress level, time and environment; 

*5. Soil may be very sensative to distur- 
bance. (Lambe, 1973b) 

It is easy to understand, therefore, the difficulties the 
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geotechnical engineer may experience in making predictions. 
Nevertheless, because of the consequences which may result 
from poor predictions, those difficulties must be overcome. 

In 1973, T.W. Lambe classified predictions accord- 
ing to the time when they are made and the information known 
at the time. This classification is shown below: 



Prediction 

Type 


When Prediction 
Made 


Results at Time 
Prediction Made 


A 


Before Event 


— 


B 


During Event 


Not Known 


B 1 


During Event 


Known 


C 


After Event 


Not Known 


C 1 


After Event 


Known 



Type A predictions are made prior to construction using only 
data available at the time. These forecasts are the challenge 
of the geotechnical engineering profession. Type B predictions 
may be based on some information obtained during the initial 
phases of construction. They are slightly less valuable than 
Type A predictions and somewhat easier to make. Class C 
predictions are "autopsies", being made after the event. 
Prediction Types B 1 and C 1 are certainly easier to make than 
the others and are of less importance. Knowing the results 
before hand may either consciously or subconsciously effect the 
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engineer's prediction procedure. Type C 1 predictions aided 
by an experienced eraser often creates a "pseudo engineering 
judgement" factor which yields amazingly accurate results. 

1.2 OBJECTIVE 

The objective of this thesis is to investigate the 
feasibility of using standardized methods to predict settle- 
ment. By making uniform the manipulation of method with para- 
meters, the author has attempted to standardize the engineering 
judgement used in each analysis. In order to evaluate these 
techniques, Type C predictions of settlement were made for 
seven case studies and then compared with actual performance 
measurements. Because Type C predictions were made, the pos- 
sibility of any "pseudo engineering judgement" was eliminated. 

A description of each one of the field cases is given in the 
next section. 

As a collateral goal, the author has tried to evalu- 
ate the role of engineering judgement in geotechnical pre- 
dictions. If engineering judgement is required in soil-rock 
predictions, at what stage should it be used and should it be 
used in conjunction with a standardized method? 

1.3 DESCRIPTION OF CASE STUDIED 
1.3.1 Alibey Test Embankment 
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The Alibey embankment was built as part of the 
Alibey Reservoir Project near Istanbul, Turkey. Initially 
built in 1967 and 1968 to evaluate performance predictions, 
the embankment was designed as a cofferdam to divert the flow 
of the Alibey Creek during the construction of the main em- 
bankment of the earth dam. The site is in a wide valley under- 
lain by alluvium consisting of about 30 meters of soft plastic 
clay with occasional sand lenses resting on a base of stiff 
clays and dense sands. The test embankment was constructed 
to a height of about 7.5 meters with the base being about 81 
meters wide and with sides sloping at about 1 on 3 . There 
are two clusters of instruments itfhich include 10 hydraulic 
piezometers, 12 vibrating wire piezometers, and 2 settlement 
plates (Dugan 1968) . The embankment and subsurface profile 
are shown in Chapter 3. 

1.3.2 Interstate Route 95 (Station 233+50) 

This section of Interstate Route 95 was constructed 
in 1968 and 1969 and is located in the town of Saugus, Massa- 
chusetts. The 35 foot high embankment is built over a low 
tidal swamp with a mean ground elevation of +5 feet and, with 
maximum high tides at about elevation +10 feet. The subsoil 
consists of a layer of peat, one of fine sand to silty sand, 
and then about 140 feet of soft to very soft gray clay which 
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is underlain by glacial till and bedrock. The gray clay 
is overconsolidated by dessication to a depth of about 80 
feet. The layer of peat directly under the 217 foot embank- 
ment base has been removed and replaced with the granular 
embankment material. The embankment sides slope at approxi- 
mately 1 on 2 (Storch Engineers, 1965). The instrumentation 
consists of 1 settlement plate, 3 hydraulic piezometers, and 4 
slope indicators. The soil profile and embankment cross 
section are shown in Chapter III. 

1.3.3 Interstate Route 95 (Station 263+00 ) 

This section of Route 1-95 is similar to that located 
at Station 233+50 with the following exceptions: the embank- 

ment rises 40 feet above ground surface and is 245 feet at the 
base; the stratum of soft to very soft clay is 90 feet thick 
(Storch Engineers, 1965). The instrumentation is similar to 
that at Station 233+50. 

1.3.4 Interstate Route 95 (Station 276+00 ) 

Again, this section of Route 1-95 is similar to 
that at Station 233+50 with the following exceptions: the 

clay layer is only about 54 feet thick, the embankment was 
constructed to about 20 feet above the ground surface; and 
it is 260 feet wide at the base (Storch Engineers, 1965). 
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The instrumentation consists of 1 settlement plate and 2 
hydraulic piezometers. 

1.3.5 M.I.T. Test Section 

This highly instrumented section of embankment is 
located at Station 246+00 of Interstate Route 95, Saugus, 
Massachusetts. The general geological features are similar 
to those described in Section 1.3.2. The embankment was 
constructed to a height of about 35 feet above the ground' 
surface with a base dimension of about 240 feet. The sides 
slope at approximately 1 on 2. The clay stratum, which 
progresses in consistency with depth from firm to soft to 
very soft, is about 136 feet thick (Guer.tin, 1962). The 
instrumentation consists of the following: 33 hydraulic 

piezometers, 6 vibrating wire piezometers, 5 observation 
wells, 12 settlement rods, 5 settlement platforms, 6 incli- 
nometers, 1 benchmark, and 9 total stress cells (Wolfskill 
and Soydemire, 1971). 

1.3.6 Northeast Test Embankment 

The Northeast Test Embankment, constructed in 1957, 
was designed as part of the Revere Beach connector to Inter- 
state Route 95 in the town of Revere, Massachusetts. The 
embankment section is approximately 60 feet long by 200 feet 
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side and rises about 35 feet above the existing ground 
surface. The subsurface conditions are similar to those 
previously described for the sections of Route 1-95 with the 
following exception: the tope 8 feet of the clay stratum is 

hard while the lower 91 feet of clay progresses in consistency 
from medium to soft. The section is instrumented with 3 settle- 
ment platforms, 44 surface displacement stakes, 4 hydraulic 
piezometers, and 2 settlement anchors (Karlsrud, 1968). 

* 

1.3.7 Lagunillas Preload 

The Lagunillas Preload project is a series of cir- 
cular embankments designed to increase the undrained shear 
strength of the foundations of three heavily loaded tanks. 
Previous investigations indicated that if the tanks were con- 
structed and filled without preloading, failure \\rould result. 

The subsoil conditions which exist at Lagunillas, on the east 
coast of Lake Maracaibo, Venezuela, consist of a stratum of 
fine silty sand which grades with depth to a sandy silt and 
finally to a soft clay. Below the clay layer lies dense 
sands and stiff clays. The particular preload studied here is 
for a tank with a diameter of 50 feet and a height of 56 feet. 
The embankment itself rises about 35 feet above ground sur- 
face, has a base diameter of 138 feet, and has sides which 
slope at about 3 on 4. Instrumentation consists of 7 
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piezometers and 5 settlement plates. The embankment and 
subsurface profile are shown in Chapter 3. 
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CHAPTER II 

METHOD OF SETTLEMENT ANALYSIS 



2.1 GENERAL 

The procedure used to make settlement predictions 
is similar to stated in Chapter I. First, the subsoil pro- 
file had to be determined as accurately as possible using 
the results of borings. Next these results were simplified 
into an accurate model of the subsurface situation. Then 
the load shape and loading time period were simplified and 
modelled. The exact standardized procedure used in modelling 
is explained in further detail in Appendix A. 

After the idealization of the field situation, the 
mechanism was determined. A mechanism is closely related to 
the method of analysis because the method, in theory, is a 
numerical model of the actual field mechanism. In each of 
the seven case studies presented here, various methods were 
used and therefore, different mechanisms were considered. 

An example of this is the Lagunillas preload case study in 
which a circular load was placed on soft clay. The situation 
is truly not one-dimensional compression, however, one method 
of analysis used was the One-Dimensional Method. Although a 
method may not be an exact theoretical representation of the 
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actual mechanism, it may yield a correct prediction when 
combined with a particular set of parameters. Additionally, 
there may be no simple method which describes the real field 
situation and an inexact method may have to be used. 

After a method was chosen, parameters were deter- 
mined from laboratory or field tests. Because the method and 
data are often incomplete and inaccurate representations of 
actual situations, and because nearly all soil engineering 
problems are statically indeterminate, methods and data 
must be manipulated to produce a reasonable prediction 
(Lambe, 1973a). The methods and paramters used to predict 
settlement for the seven cases studied are described briefly 
in the following sections. More detailed, step-by-step 
techniques are listed in Appendix B. 

2.2 METHODS AND PARAMETERS USED TO PREDICTE TOTAL SETTLEMENT 
2.2.1 One-Dimensional Method 

This method assumes that the mechanism which causes 
settlement is one-dimensional strain. It also assumes that 
the change in vertical stress is constant with depth and that 
no undrained deformation (initial settlement) occurs within 
the cohesive strata. Initial settlement due to loose cohe- 
sionless layers is computed, however. The instances in which 
one-dimensional strain occurs may be rare. Nevertheless, 
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when the loaded area becomes very large and the thickness 
of the compressible layers small, the strains under the 
center of the load are then closely approximated by this 
mechanism. 

The parameters required for this type of analysis 
are obtained from the void ratio versus logarithm stress 
(e-log p) curves of oedometer tests on samples from various 
depths throughout the cohesive layers. Schmertmann ' s method 
for overcoming the effects of sample distrubance should be 
used to determine the maximum past pressure (o" vm ) and to 
construct the "field" e-log p curve (Ladd, 1971). The values 
for the Compression Ration (CR) and the Recompression Ratio 
(RR) are taken from the field curve and a weighted average of 
each is found over the depth. The initial effective vertical 
stress (cT ) should be determined by the procedure set forth 
by Lambe and Whitman (1969) . The final vertical effective 
stress plus the applied load. 

The One -Dimensional Method requires that compressible 
strata be divided into layers and settlements be computed for 
each layer. Recommendations for layering of the subsoil are 
located in Appendix A and the computational techniques in 
Appendix B. 
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2.2.2 Modified One-Dimensional Method 

This method is very similar to the one pre- 
viously described. The only major difference being that 
the change in vertical stress (Ag v ) is not constant with 
depth and thus the final vertical effective stress (cf^) 
differs from that of the One-Dimensional Method. The 

* 

change in vertical effective stress with depth should 
be computed from elastic theory where Poisson's Ratio is 
equal to one-half (y^=l/2) . Charts of elastic solutions 
have been published by Lysmer and Duncan (1972) , Davis 
and Poulos (1973), Lambe and Whitman (1969), and others. 

The final effective vertical stress is the sum of the 
initial vertical effective stress and the change in vertical 



stress . 
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2.2.3 Elastic Method 

The Elastic Method assumes that a soil stratum is 
an isotropic elastic-medium which can be represented by a 
set of elastic constants. When the constants have been deter- 
mined, predictions can be made with the aid of charts and 
tables which have been published (Lysmer and Duncan, 1972; 

Davis and Poulos, 1968). Total settlement is considered to 
be composed of two components, initial elastic settlement 
(p ) and consolidation settlement (p ), due to undrained ^nd 
drained conditions resprectively . Both components are com- 
puted by the same method but different paramters. The para- 
meters required for initial elastic settlement predictions 
are an undrained Young's Modulus (E^) and an undrained Poisson's 
Ration (y^). The parameters required for computation of con- 
solidation settlement are a drained modulus (E ) and a drained 
value of Poisson's Ratio (y c ) . The procedure for determina- 
tion of these parameters and preferences recommended in 
selecting the data are given in Appendix B. When values of 
Young's modulus vary with depth, a weighted average should 
be determined for the analysis. 

The techniques used in computation of the Elastic 
Method are taken from a compilation of tables and charts of 
elastic solutions published by Lysmer and Duncan (1972). 

Other sources are available, but all Elastic Method 
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predictions made for the seven cases studied here were 
computed by means of the Lysmer and Duncan solutions. 

2.2.4 Elastic Method Considering First Local Yield 

This method, proposed by D'Appolonia, Poulos, 
and Ladd (1971) , modifies the inital settlement predicted 
in Section 2.2.3 to account for yielding of highly stressed 
areas. When first local yield occurs, the stress - strain 
properties of that volume of soil become nonlinear and plastic 
flow occurs. As a result, the initial undrained deformation 
is no longer linearly elastic. This method utilizes graphs 
prepared from the results of nonlinear finite element analyses 
to give a factor, the Settlement Ratio (S r ) , which can be 
applied to modify the initial eastic settlement (P ) 
(D'Appolonia, et al . , 1971). In order to use these graphs, 
the following parameters must be known: 

1. The loading and subsurface 
geometry (H/B) ; 

2. The ultimate bearing capacity 
(q u i t ) which is a function of 

. the undrained strength (s u ) and 
geometry; 

3. The intial shear stress ration ( f) 
which is a function of the initial 
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vertical and horizontal effective 
stresses (cT V0 and ^ ) and the un- 
drained shear strength; 

4. The inital elastic settlement C P e ) 
computed by the method described in 
Section 2.2.3. 

2.2.5 Skempton-Ej errum Method 

The basic premise of this method is that consolida- 
tion settlements occur as a result of the dissipation of 
excess pore pressures (Au) and that the excess pore pressures 
are not necessarily equal to the change in vertical stress 
caused by the applied load. Additionally, the Skempton- 
Bj errum Method allows for initial settlements due to undrained, 
non-one-dimensional , deformation. Initial settlements are 
computed from either the purely elastic method (Section 2.2.3) 
or the method which considers the effect of first local yield 
(Section 2.2.4). 

Since consolidation settlement is a function of the 
excess pore pressure, an important paramter required for this 
method is Skempton's A parameter. Also, because the dissi- 
pation of the pore pressure is what causes the settlement, 
the loading geometry and stratigraphy must be considered. 
Skempton and Bjerrum (1957) proposed a semi- empirical factor 



